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ABSTRACT 

Determining magnetic field properties in different environments of the cosmic large-scale structure as well as their evolution over 
redshift is a fundamental step toward uncovering the origin of cosmic magnetic fields. Radio observations permit the study of extra- 
galactic magnetic fields via measurements of the Faraday depth of extragalactic radio sources. Our aim is to investigate how much 
different extragalactic environments contribute to the Faraday depth variance of these sources. We develop a Bayesian algorithm to 
distinguish statistically Faraday depth variance contributions intrinsic to the source from those due to the medium between the source 
and the observer. In our algorithm the Galactic foreground and the measurement noise are taken into account as the uncertainty corre¬ 
lations of the galactic model. Additionally, our algorithm allows for the investigation of possible redshift evolution of the extragalactic 
contribution. This work presents the derivation of the algorithm and tests performed on mock observations. With cosmic magnetism 
being one of the key science projects of the new generation of radio interferometers we have made predictions for the algorithm’s 
performance on data from the next generation of radio interferometers. Applications to real data are left for future work. 

Key words. Methods: data analysis - methods: statistical - magnetic fields - polarization - large-scale structure of Universe 


1. Introduction 

The origin and evolution of cosmic magnetism are at present 
poorly understood. Answering the many open questions sur¬ 
rounding the physics of astrophysical magnetic fields is a dif¬ 
ficult task since magnetic fields can be significantly affected by 
structure and galaxy formation and evolution processes. Their 
strength can be amplified for example in galaxy clusters by 
mergers, and in galaxies by large-scale dynamos, invoking differ¬ 
ential rotation and turbulence. Insights into the origin and prop¬ 
erties of magnetic fields in the Universe could be provided by 
probing them on even larger scales. Along filaments and voids 
of the cosmic web, turbulent intracluster gas motions have not 
yet enhanced the magnetic field; its strength thus still depends 


on the seed field intensity, in contrast to galaxy clusters, where it 
probably mostly reflects the present level of turbulence (see e.g., 
ponnert et al.|2009l|Xu et aL|2010||201 l| l. Intervening magneto¬ 
ionic media cause a difference in fhe phase velocity between the 
left-handed and right-handed circular polarization components 
of the linearly polarized synchrotron radiation emitted by a back¬ 
ground radio source (e.g. |Carilli & Taylor|2002[ |Govoni & Fer-| 
|etti ||2004] l. This effect translates into a rotation of the intrinsic 
polarization angle, if/Q, 

^ il,o + ( 1 ) 


Following |Burn| ( 'l966| l, the observed polarization angle, i/r, de¬ 
pends on the observation wavelength, A, through the Faraday 
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depth, (p. 


(p^ao 


r 


dl 

Bi(z) n^iz) dz, 
dz 


( 2 ) 


where gq depends only on fundamental constants, is the elec¬ 
tron density, Bi the component of the magnetic field along the 
line of sight, and Zs the redshift of the source. When the rotation 
is completely due to a foreground screen, the Faraday depth has 
the same value as the rotation measure (RM), defined by 


RM = 


dtp 


( 3 ) 


The Faraday depth is assumed to be positive when the line of 
sight average component of the magnetic field points toward the 
observer, otherwise it is negative for a field with an average com¬ 
ponent pointing away from the observer. The amount of Faraday 
depth measured by radio observations along a given line of sight 
is the sum of the contributions from the Milky Way, the emit¬ 
ting radio source and any other sources and large-scale struc¬ 
tures in between hosting a magnetized plasma. The investiga¬ 
tion of these contribution and of their possible dependence on 
redshift is essential in order to discriminate among the different 
scenarios of magnetic field formation and evolution and there¬ 
fore crucial for the understanding of cosmic magnetism. Sen¬ 
sitive observations, a good knowledge of the Galactic Faraday 
foreground screen, and a statistical approach able to properly 
combine all the observational information are necessary. An all¬ 
sky map of the Galactic Faraday rotation foreground and an es¬ 
timate of the overall extragalactic contribution has been derived 
by |Oppermann et ^(|2012||2015[) in the framework of Informa¬ 
tion Field Theory ( EnBlin et al.|2009 1 , by assuming a correlated 
Galactic foreground and a completely uncorrelated extragalac¬ 
tic term. In this paper, we propose a new, fully Bayesian, ap¬ 
proach aiming at further disentangling the contribution intrinsic 
to emitting sources from the contribution due to the intergalactic 
environment between the source and the observer, and at inves¬ 
tigating the dependence of these contributions on redshift. 

The first direct proof of the existence of magnetic fields 
in large-scale extragalactic environments, i.e., galaxy clusters, 
dates back to the 1970s with the discovery of extended, diffuse, 
central synchrotron sources called radio halos (see e.g. Feretti et 
|al.|2012| for a review). Later, indirect evidence of the existence of 
intracluster magnetic fields has been given by several statistical 
studies on the effect of the Faraday rotation on the radio sig¬ 
nal from background galaxies or galaxies embedded in galaxy 
clusters ( [Lawler & Dennison|1982[|Vallee et'd? 1986[|Cl^ke et 


ay|200]J Johnston-Hollitt||2003| |Clarke||2004[ [johnston-Hollitt 

& Lkers||2004j l. On scales up to a few Mpc from the nearest 
galaxy cluster, possibly along filaments, only a few diffuse syn¬ 
chrotron sources have been reported (Harris et al.|1993 [Bagchi 


|et al.|20d2 [Kronberg et al.|20W)|Giovannini et ai|2013[|2015| l. 

Magnetic fields with strengths of the order of 10^^^ G in voids 

might be indicated by y-ray observations (see [Neronov & Vovk 


2010[|Thvecchio et al.|2010l|Takahashi et al.|26l2||2013| but see 

Broderick et al.|2014a|b for alternative possibilities). Neverthe¬ 
less, up to now, a robust confirmed detection of magnetic fields 
on scales much larger than clusters is not available. [Stasyszyn] 
jet al.| ( |2010| l and jAkahori et al.] ( |2014a| ) investigated the possi¬ 
bility to statistically measure Laraday rotation from intergalac¬ 
tic magnetic fields with present observations, showing that only 
the Square Kilometre Array (SKA) and its pathfinders are likely 
to succeed in this respect. By comparing the observations with 
single-scale magnetic field simulations, jPshirkov et al.j P015| l 
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infer an upper limit of 1.2 nG for extragalactic large-scale mag¬ 
netic fields, while the [Planck Collaboration j ( [2015| l derived a 
more stringent upper limit for primordial large-scale magnetic 
fields of Z? < 0.67 nG from the analysis of the CMB power spec¬ 


tra and the effect on the ionization history (but see also Takahashi 
[et al.|2005[[Ichiki et al.|2006t . 

A number of authors examined a possible dependence of ex¬ 
tragalactic Laraday depths on the redshift of the observed radio 


source, but no firm conclusion has yet been drawn. Kronberg & 


[Perry[ ( [T982 | i found an increased variance of the Laraday depth 
in conjunction with higher redshifts, as also found in some later 
studies (e.g. Welter et al.|1984[[Kronberg et al.|2008 1 . However, 
[Oren & Wolfe ( T995| l did not find any evidence of an increase of 
the variance of the Laraday depth as a function of the redshift, 
as suggested also by the recent work by Hammond et al. ( 2012| l 
and [Pshirkov et ari ( [2015[ l. 

The rest of this paper is organized as follows. In §[^ we 
describe the theory behind our method. In §[^ the tests per¬ 
formed are presented, with predictions for the new generation 
of radio interferometers. Moreover, we outline a generaliza¬ 
tion of the algorithm in order to discriminate the contribution 
from different large-scale structures along the line of sight. Li- 
nally, in §[^ we draw our conclusions. The application of the 
algorithm to real data is left for future work, as explained in 
the text. In the following we adopt a ACDM cosmology with 
Hq - 67.3kms 'Mpc ', Qni=0.315, Qa= 0.685, and Qc=0.0 
( [Planck Collaboration 12014] ). 

2. Theoretical framework 

The probability for the extragalactic contribution, along the 
line of sight, i, to take on a value within the infinitesimal inter¬ 
val between 0e,i and cp^j + given the data, d, and a model 
parameterization, 0, is given by the probability density distri¬ 
bution. In turbulent environments (see, e.g., Lawler & Dennison] 


[T982l|Tnbble| 1991 1 [Feretti et al.|1995||Lelten|1996[ for galaxy 

clusters), this probability density distribution can be represented 
with a Gaussian. The mean and variance. 


of this distribution are 


= 0 
((Pli) = 4 


ff 


( 4 ) 


( 5 ) 


^ ^dzdz'(n(z)n(z')B,(z)Bi(z'). 
0 dzdz' 


The zero mean results from the fact that we do not have any rea¬ 
son to suppose either a positive or a negative mean Faraday depth 
value. Frozen-in magnetic fields are expected to have strengths 
depending on redshift z. Because of the expansion of the Uni¬ 
verse, lengths are stretched 

lo 


= (1+Z)- 


( 6 ) 


In an isotropically expanding fluid with no significant fluctua¬ 
tions, the mean electron density evolves as 


(no) 


( 7 ) 


and, if the magnetic flux density is conserved, the magnetic field 
strength as 


./<!!£> V . n + 
(Bo) \(no)j 


( 8 ) 
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Fig. 1. Distribution of the real and mock sample of redshifts z {top), 
and flux densities Fjy at 1.4 GHz {middle). In the bottom panel, flux 
density versus redshift for the sources with both measurements available 
( [Taylor et al.|2009[|Hammond et al.|2012[ l. 


where Iq, (no), and {Bq) are the present-day values, and /, (ng), 
and {B) are the values at the time when the signal was emitted 
by the source. If we consider these assumptions and we define a 


length scale A/ := f dl'{B{l)B{l'))I{B{r)^), we obtain 

<<,) - «o j\nl){Bl}A,o{l + (9) 

For the derivation of this expression see Appendix[A| Here, we 
assumed an unstructured Universe and used the definition of 
proper displacement along a light-ray, ^ = (i+;)Hfe) ■ Moreover, 
we assumed that, within a correlation length A;o, the redshift 
can be approximated to be constant. Magnetic field strength and 
structure as well as the electron density have different values in 
different environments, j. In the following, we assume them to 
depend only on the environment and not on the location within 
an environment. This simplification renders the problem feasi¬ 
ble. 

In this paper, as a first step, we will restrict our analysis to a 
two components (scenario 2C), the emitting radio source itself, 
whose contribution is o-^^^.{Zi,&), and the medium between the 
source and the Galaxy, whose contribution is -{zi, ©), such 
that 


0) = ®) + 0)- 


( 10 ) 


We denote with © = {crint,o, CTenv.o^Uium.Ared) an A-dimensional 
vector, where N is the number of parameters used in the repre¬ 
sentation of the Faraday depth variance. These are the parame¬ 
ters we want to infer. We choose the following parameterization 
for the intrinsic contribution to the variance in Faraday depth, 

.^2 




Li 


int.O 


. ( 11 ) 

(l+Z/)4 

where L,- is the luminosity ^ of the source i, Lq = 10^’W/Hz, 
and xium absorbs possible dependencies on the luminosity of the 
source since faint sources may not be detected, and for the envi¬ 
ronmental contribution 




L^iiZi,Xred} 


where Dq - 1 Gpc, and Di{z,Xied) is defined as 


f 


-( 1 - 1 - z)'^^^“'‘dz, 


( 12 ) 


(13) 


DiiZi,Xred) ■ . 

Jo H{z) 

to capture the redshift scaling of Eq. (|^ as well as auxiliary mod¬ 
ifications by the dumpiness of the Universe Vm Xmd- Since the 
length of the path covered by the signal in the source is unknown, 
we factored it in g in Eq. 0 We note that (Xint o and cr^av.o 
have been assumed to be independent of the redshift. In Eq. ( fTT) 
the only redshift dependence is absorbed by the factor (1 -H ZiT^ 
that takes into account the effect of redshift on Earaday rotation 
(squared), while in Eq. dl^ it is absorbed by Di{zi,Xied) which 
takes into account Eq. (|13[). 

This parameterization describes the simplest scenario. Eor 
more complex scenarios that include three components, we re¬ 
fer to Appendix^ where we introduce an additional constant 
and latitude dependent term. Moreover, galaxies along the line of 
sight between a source and the observer can be responsible for 
high Earaday depth values (e.g. Kronberg & Perry] 1982[ [Wel-| 


ter et al.||l984 Kronberg et al.|[^08t Bernet et al.||2012| l, indi¬ 

cating magnetic field strengths in these intervening galaxies of 
(1.8+0.4) juG ( [Earnes et al.|2014^ . The rotation of the polariza¬ 
tion angle due to these sources adds to that associated with the 
large-scale structure and, therefore, should be taken into account 
in a proper modeling. Since the aim of this paper is to give a 
proof of concept, we leave this for future work. 


* This luminosity refers to the mean-frequency of the frequency band 
used for the computation of the Faraday depth values. 
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2.1. Bayesian inference 

To constrain the vector 0 = {cfint.o. o-env.o^Tium.Tred) on the basis 
of these data, d, we propose a Bayesian approach. The posterior 
probability distribution, P{s\d), on a signal, s, after a dataset, d, 
is acquired, can be expressed with Bayes’ theorem. 


P(s\d) = 


P{d\s)P(s) 

P(d) 


(14) 


The prior probability distribution, P(s), is modified by the data, 
d, through the likelihood, P{d\s). The evidence, P(d), is a nor¬ 
malization factor, obtained by marginalizing the joint probabil¬ 
ity, P(d, s) = P(d\s)P(s), over all possible configurations of the 
signal, s. 

In this context, the data, d, can be represented as a vector 
with elements di, with 1=1,..., A^ios, where Nios is the total num¬ 
ber of lines of sight. Each measurement, dj, is the Faraday depth 
evaluated in the direction of the source, i, and is the result of the 
sum of a Galactic and an extragalactic contribution, 0g_, and ^e,/, 
and the noise, n,, of the measurement process,such that 


di - 


■m. 


(15) 


From the observed data, the Galactic foreground should be re¬ 
moved as well as possible to reveal the extragalactic contri¬ 
bution. However, any estimation of the Galactic foreground is 
based on the same data and is facing the separation problem for 
Galactic and extragalactic contributions. The only available dis¬ 
criminant (so far) is the large angular correlation the Galactic 
contribution shows. This allow for a statistical separation and 
Galactic model construction. Such a model will inevitably have 
uncertainties and correlations among such uncertainties, which 
have to be properly taken into account in a statistical search 
for extragalactic Faraday signals. To this end, |Oppermann et al.| 
( 2012| 2015| l developed a fully Bayesian approach to reconstruct 
the Galactic Faraday depth foreground and estimate the extra¬ 
galactic contribution as well as the involved uncertainties by us¬ 
ing the Faraday depth catalogs available to date. Their posterior 
for the extragalactic contribution can be used in order to further 
disentangle the intrinsic and environmental contributions. |Op-| 
permann et al. (2012 2015| l’s analysis relies on the assumption 
that for each source the prior knowledge can be described by a 
Gaussian probability density distribution. 


Nlos 

Pirp.) (16) 

with a standard deviation, erg =» 7 radm^^, irrespective of the line 
of sigh0 Here, on the other hand, we want to test whether the 
variance is different for each line of sight, depending on the red- 
shift of the source, according to Eq. 0 and Eq. ( [T^ . We note 
that in our approach the angular separation of sources is assumed 
to be large enough that the magnetic fields probed by different 
lines of sight can be modeled as uncorrelated. All angular corre¬ 
lations of Faraday depth on scales down to the effective resolu¬ 
tion of the catalog (a!0.5°), are absorbed by the Galactic compo¬ 
nent in this model. A more complete modeling would require to 
take into account possible correlations in the extragalactic com¬ 
ponent. Thus, our assumptions imply that the values derived with 
the proposed algorithm for the contributions of different environ¬ 
ments to the Faraday depth dispersion are a lower limit. 

^ The notation Q(x, X) indicates a one-dimensional Gaussian distribu¬ 
tion for a variable x with zero mean and variance X. 
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We cannot follow the prescription described in Appendix 
D.2.3 of [Oppermann et ’ak] ( 2015] l because our prior assump¬ 
tions are too different from theirs. Instead, we resort to Gibbs 
sampling ( |Geman & Geman~||1984[ [Wandelt et al.|[2004) l. This 
approach relies on the fact that sampling from the conditional 
probability densities. 




(17) 


and. 


0^B(0|<^e,t!?), (18) 

in a two-step iterative process is equivalent to drawing samples 
from the joint probability density 

^ P(®,<P^d), (19) 

if the process is ergodic. 

For the parameters cTiatfl and crgnv.o we choose a prior, 

P{o-mifl, cre„^_o) const. (20) 

Results with other priors are discussed in Appendix[D] Con¬ 
versely, we do not expect the parameters xium and Xred to differ 
greatly from 0, since we have already accounted for all obvious 
redshift effects. This requirement is satisfied if we use the fol¬ 
lowing Gaussian priors 


PiXlnm, Tred) = &(Xlum, 1 WOfred , 1), (21) 

and in combination 


P(0) = P(crint,o, Cr env, o)P(xi um? Xred)- 


( 22 ) 


2.2. Description of the aigorithm 

Here, we describe the Gibbs sampling procedure mentioned in 
the previous section. We run the algorithm starting from values 
of the parameters 0 randomly drawn from their prior. This 0 
vector is used to compute a variance for the prior of the extra¬ 
galactic contribution. 


N,„, 

Pirp.m = f]^(0e,,-,fre.,(2-0))> 


(23) 


1=1 


where cr^ ^.(z;, 0) is evaluated according to Eq. (10 1 . A sample for 
the extragalactic contribution is drawn from the posterior 


P{^.\®,d) 


p{d\cp„®)p{^,m 

P{d\Q) 


(24) 


following the approach described in [Oppermann et ar] ( |2015) l, 
with the Galactic power spectrum, the Galactic profile, and the 
correction factors to the observed noise variance (indicated by rp 
in the paper by [Oppermann et al.||2015| l fixed to the published 
values. After fixing the extragalactic sample, a new 0 sample is 
drawn from the conditional probability 


P{&\4>.) = 


R(0e|Q)R(Q) 

P(A) 


N,„. 

1=1 


(25) 


Here, we drop the dependence on the data, d, because 0 and 
d are conditionally independent given (j)^. To sample from this 
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Table 1. Mock catalogs. Col. 1: scenario considered for the generation of the mock catalog. Col. 2: identification code (ID) of the test. Col. 3: 
short description of the catalog. Col. 4: number of lines contained in the catalog. Col. 5: Observational uncertainty in Faraday depth. Col. 6: Figure 
where the results of the corresponding test are shown. 


Scenario 

ID 


Description 




los # 

cr^ . 

noise 

(rad/m^) 

Figure 

Two Component 

2C1 

Datasets used by 

Oppermann et al. 

(2012 

2015 


41632 

13.0 


^a) 



2C2 

Datasets used by 

Oppermann et al. 

(2012 

2015 

1* 

4003 

13.0 


^b) 



NPC 

LOEAR 120 

-160 MHz, North 

“"olar Cap 


2148 

0.05 


a,b) 


GW 

LOFAR 120-160 MHz, Great Wall 


1036 

0.05 


E.l 



POSSUM 

ASKAP-POSSUM 1130-1430 MHz, South Polar Cap 

3476 

6.0 





B2SPC1 

SKA Band 2, 0.65-1.67 GHz, South Polar Cap 

3476 

0.8 

4 


) 


B2SPC2 

SKA Band 2, 0.65-1.67 GHz, South Polar Cap 

1129 

0.8 


I'h) 


BlSPCl 

SKA Band 1, 0.35-1.05 GHz, South Polar Cap 

3476 

0.3 

K 

r a,b) 


B1SPC2 

SKA Band 1, 0.35-1.05 GHz, South Polar Cap 

1129 

0.3 

F 

z a) 

Three Component 

3C1 

Datasets used by 

Oppermann et al. 

(2012 

2015 

1* 

41632 

13.0 

C 

TTa) 


3C2 

Datasets used by 

Oppermann et al. 

(2012 

2015 


4003 

13.0 

C 

.1 b) 

Latitude Dependence 

EDI 

Datasets used by 

Oppermann et al. 

(2012 

2015 


41632 

13.0 

C 

^ ’a) 


LD2 

Datasets used by 

Oppermann et al. 

(2012 

2015 


4003 

13.0 

C 

:2 b) 

Prior 0 

PO 

Datasets used by 

Oppermann et al. 

(2012 

2015 


41632 

13.0 

D 

TTa) 

Prior 1 

PI 

Datasets used by 

Oppermann et al. 

(2012 

2015 


41632 

13.0 

D 

.lb) 


“ For present observations this value represents the mean value of the observed uncertainties.The references for the surveys and catalogs used 


by Oppermann et al. ( 2012^ are: Dennison \\919 1; 

labara & Inoue 1 1980|l; Simard-Normandin et al. 

l|198ip; Lawler & Dennison|) 1982|l; Rudnick 

& Jones|11983 i;|Kato et al. (1987 

|;|Broten et al.|i 

1988 i;|Hennessy et al.|11989J; Kim et al. ( 1991J; 

Clegg et al. 1 1992 1 ; Wrobel|11993 i; Oren & 

Wolfe|(|1995);|Minter &. Spangler 

|1996J; Condon et al. 11998J; Gregorini et al. ( 1998J; 

Vigotti et a 


Clarke et al.|( 2001|(;|Gaensler et al. 

(|2001J;|Brown et al. 

1 2003 i; Johnston-Hollitt 120031;|Klein et al. 

12003 

i;|Taylor et al. 

20031; (^^lar 

lcelj2004i 

Johnston-Hollitt & Ekers|12004 1 ; 

Gaensler et al. 1 2005 

; McClure-Griffiths et al.|l2005|l; Roy et al. 

2005 

Haverkorn et al. 12006J; [Brown et a 

l 1.|( 2007|I; Braun et al.|12007 1 ; Mao 

et al. 12008 1 ; Feain et al.|( 2009 1 ; Heald et al. 1 2009f; Taylor et al. 

2009 

; Bonafede et al. 12010J;|Mao et al.|( 2010 1 ; Feain et al. ( 201H;|Van Eck 

et al.|i 

mTIT 


distribution we use a Metropolis-Hastings algorithm ( |Metropo^ 
lis et al. ||1953[ Hasting || 1970 1. When direct sampling is diffi¬ 


cult, Metropolis-Hastings algorithms can approximate a proba¬ 
bility distribution with random samples generated from the dis¬ 
tribution itself. At each iteration a step in the parameter space 
is proposed according to a transition kernel and then accepted 
according to an acceptance function. If the proposed step is not 
accepted, the old 0 values are kept and used to draw a new sam¬ 
ple of the extragalactic Faraday depths. 

The convergence criteria adopted in this work are described 
in Appendix [B| 


3. Results 


In the following we present tests performed with different Fara¬ 
day depth catalogs. These catalogs differ in the number of com¬ 
ponents used to generate the overall extragalactic Faraday depth, 
the numbers of lines of sight in the sky, and the observational un¬ 
certainties that are different for the different radio surveys con¬ 
sidered here. In Sect. we demonstrate that the algorithm is 
working properly for the two-component scenario described in 
Sect.|^ In Sect. |3.2| we present the prospects with the surveys 
planned with the new generation of radio interferometers (LO- 
FAR, ASKAP, and SKA). 

We perform tests assuming an overall extragalactic Faraday 
depth in agree ment with the values presently inferred, na mely 
ss 7.0rad m“^ ( S chnitzel^|20 10 Oppermann et al.||2015 i, and 
comparable intrinsic and the environmental contributions. To 
satisfy these two requirements, we need to use slightly different 
values of the © parameters for surveys with different frequency 
specifications. Indeed, the contributions depend on the frequency 
through the luminosity of the source, see Eq. O- Our choice of 


the © parameters translates to a strength of magnetic fields in¬ 
trinsic to the source of 


(Bio) 

liG 


~ 0.5 H- 1 


(«o) ^ ‘ / A;o \ V L \ 

10 “^ cm~3 / \ 5 kpc / \ 100 kpc / 


(26) 


where L is the size of the emitting radio source, and to large scale 
magnetic field strengths of 


{Bm) ^1 <«o) 'i V A/o \ 

nG \10“^cm“3/ \5Mpc/ 


(27) 


In the tests for LOFAR, ASKAP, and SKA, we additionally con¬ 
sider an overall extragalactic Faraday rotation of ^ 0.7 rad m“^ to 
mimic weaker fields. This corresponds to magnetic field values 
weaker by a factor ten than those given in Eq. ( p^ and Eq. ( |27] l. 

The two-component parameterization represents the simplest 
scenario. Nevertheless, the algorithm is able to successfully deal 
also with more complex scenarios that include a third constant 
or latitude-dependent component. These scenarios are discussed 
in Appendix[C| Moreover, to asses if different priors can have an 
impact on our results, we consider in Appendix[D] a flat prior in 
cr^ and a flat prior in In(cr^). In Table[^we give a summary of all 
the setups we used, including those presented in the appendices. 
To each of them we assigned an identification code (ID), used in 
the paper to discriminate among the different scenarios. When 
the same collection of sources is used for tests with different 
values of the overall extragalactic Earaday depth, we distinguish 
among them by adding a roman letter. For example, Xa and Xb 
indicate tests performed using the collection of sources X and 
two different overall extragalactic Faraday depth standard devi¬ 
ation, a denotes 7.0radm“^, while b ^ 0.7radm^^. In Table|^ 
we report all the values of the © parameters adopted in the differ¬ 
ent tests. In the main text we give both a quantitative and a visual 
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summary of the results of all the tests we performed, while the 
full posteriors are shown only for the most important tests we 
can'ied out. To make the reading easier, the posteriors for all the 
other tests are shown in the appendices. 


3.1. Present-instrument observations 


To assess the quality of the algorithm we generate a mock cata 
log for the sample of sources used by |Oppermann et ^ (2012 


2015|l. This mock catalog includes coordinates, redshifts, lumi¬ 


nosities, and Faraday depth values. 

The positions of the sources on the sky were kept the same 
as for the real sources. The majority (a!40000 sources) belong to 
the catalog of Taylor et al.| ( |2009| l, and for 4003 of them, spec¬ 


troscopic redshift measurements have been published by Ham- 


mond et al. 


(20121. For most of the sources, these catalogs give 


a flux density measurement that allows the computation of the 
luminosity of the source. Where available, we use the measured 
redshift and flux density. For the vast majority of the sources we 
generate a mock redshift, and for a few of them a mock flux den¬ 
sity value. Mock redshifts and flux densities are extracted inde¬ 
pendently from the two observed distributions. In Fig.[2 the dis¬ 
tribution of both the real and mock sample of redshifts and flux 
densities is shown respectively in the top and middle panel. In the 
bottom panel the observed flux density versus redshift distribu¬ 
tion is presented. These two quantities appear to be weakly cor¬ 
related. For sake of simplicity we neglected such correlation in 
our mock simulation, since it should not have any impact on our 
analysis. We assume all redshifts and luminosities to be known 
with negligible uncertainty. 

For all the sources in the catalog we generate a mock Fara¬ 
day depth value. The observed Faraday depth values consist of a 
Galactic, an extragalactic, and a noise contribution. We consid¬ 
ered the Galactic contribution to be given by a sample extracted 
from the posterior of |Oppermann et JI] ( |2015 i. To mimic obser¬ 
vational uncertainties, the noise variance has been computed for 
each source according to Eq. (37) in |Oppermann et ak ( 2015| l, 
where as observed uncertainty, cr,, we use the uncertainties re¬ 
ported in the observational catalogs and as 77 , we use the val¬ 
ues recovered by |Oppermann et ak] P015| l. The observational 
error of each measurement has been extracted from a Gaussian 
with this standard deviation and zero-mean. Concerning the ex¬ 
tragalactic contribution, in this test, we consider the 2C scenario 
described in §[^ namely an intrinsic and an environmental con¬ 
tribution. 


mliizi,®) 


DiiZi) 2 

io/ (l+z,)4^ £>o 


(28) 


The variances in Eq. ( [28| depend on the redshifts and lumi¬ 
nosities of the sources. Therefore, each source will have a differ¬ 
ent variance. Eor each source, the extragalactic contribution is 
extracted from a Gaussian with this variance and zero mean. As 
summarized in Table the contributions intrinsic to the source 
and to the medium between the source and the observer are re¬ 
spectively, To-ijio = 18.2radm“^ and = 1.4radm“^. Since 
the mean of the factor L, 7 Lo/(l + 2 ;)^ is ~ 0.06 and the mean of 
the factor Z),/77o is ~ 15.5, the standard deviation in the overall 
intrinsic and environmental components are crl™ ss 4.4 rad m“^ 
and cr‘™ ai 5.3 radm“^. Eor this scenario, we run two tests cor¬ 
responding to a different number of lines of sight: 

- 41632 (2C1). This is the total number of lines of sight for 
which an estimate of the extragalactic contribution is avail¬ 
able from|Oppermann et al.|P015]l; 


- 4003 (2C2). This number accounts for all the sources in the 
catalog 2 C 1 for which a redshift measurement is available as 
well ( [Hammond et al.|2012| l. 

In Fig.[^ we show the results for the tests 2C1 and 2C2, 
meaning with a two-component scenario for 41632 lines of sight 
(a) and 4003 lines of sight (b). The histograms in the top plots of 
each column represent the 1 -dimensional projection of the pos¬ 
terior for the corresponding parameter. The dotted lines mark 
the true value of the parameter. The dashed and dashed-dotted 
lines describe the posterior statistics, namely the mean and the 
Icr confidence level respectively. The continuous lines indicate 
the prior used in our analysis. The panels in colors show the 2- 
dimensional projection of the posterior for a given couple of pa¬ 
rameters. These plots show that our algorithm is able to recover 
the mock ©-values for this scenario. The inferred posterior mean 
values agree well within the uncertainties with the correct ones 
and the posterior distributions are much narrower than the prior 
distributions. The dispersion in the parameters 0 increases by 
decreasing the number of lines of sight, as expected. The plots 
indicate that some of the parameters are correlated, e.g., most no¬ 
ticeably cTgny oWred that show a strong anticorrelation. This fea¬ 
ture can be understood in light of Eq. ( [T^ . Indeed, for a given 
Earaday rotation cr^nv associated with the structures between the 
source and the observer, larger crg„y Q imply smaller;t'red vice 
versa. We expect the correlation in the posterior to be signifi¬ 
cant for any reasonable parameterization that allows for the same 
number of degrees of freedom. 

In order to have a compact and complete visualization, in the 
rest of the paper we present all the tests we performed and their 
results as in Eig.[^ 


3 . 2 . Future prospects 


We investigate the possibility to separate the Earaday rotation 
intrinsic to the emitting radio source from that due to the ex¬ 
tragalactic environments between the source and the observer 
with the specifications of the SKA, its precursor, ASKAP, and 
its pathfinder, LOEAR. The mock rotation measure values for 
each source have been generated as described in Sect. ED Eor 
each catalog the noise contribution has been extracted for each 
line of sight from a Gaussian with zero-mean and standard devi¬ 
ation equal to maximum uncertainty in Earaday depth expected 


in the corresponding frequency range, according to Stepanov et 
( |2008| ). Luminosities at frequencies different from 1.4 GHz 
have been spectrally adjustecj^ 


Since our approach assumes that all the lines of sight are 
independent, i.e. sufficiently separated (> 1 °) so that the cross¬ 
correlation function of their magnetic field is zero, in the follow¬ 
ing tests we compute the number of lines of sight considering a 
density of sources lower than or equal to one polarized source 
per square degree. The number of sources detected by ASKAP 
and SKA per square degree will be at least one hundred times 
larger. Here, we investigate if already with a small sample of 
lines of sight we would be able to put any constraints on the 
contribution from extragalactic large-scale environments, rather 
than exploit the information delivered by the full number of lines 
of sight and to show the full potential of ASKAP and SKA. 


^ For convenience, we use a single power-law F(v) oc v where F(v) 
is the flux density. 
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3.2.1. LOFAR 


We consider the LOFAR’s High Band Antennas (HBA) only, 
because of their better polarization performance. We select the 
HBA region of the spectrum with less contamination from ra¬ 
dio frequency interference (120-160 MHz, Offringa et al. 2013| l. 
In this frequency range, the uncertainty on rotation measure val¬ 
ues is expected to be < 0.05 rad m“^ for a signal-to-noise ratio 
larger than 5 (Stepanov et al. [2008)1. According to Mulcahy et 


[ar| ( |20T^ , the expected number of polarized extragalactic ra¬ 
dio sources is 1 per 1.7 square degrees for 8 h-long observations, 
assuming an average degree of polarization of 1%, a spatial res¬ 
olution of 20", and a detection threshold of 500//Jy/beam/rms^ 
that corresponds to S/N = 5. 

On the basis of these assumptions we generate coordinates 
in the sky for a collection of sources corresponding to a sur¬ 
vey (8 h per pointing) in the direction of the North Polar Cap 
(NPC). Among these sources we select those with Galactic lati¬ 
tude larger than 55°. This results in approximately 2200 sources. 
We derive a catalog assuming an overall extragalactic Faraday 
rotation o-g xil rad m“^ (NPCa), and one assuming an overall ex¬ 
tragalactic Faraday rotation o-g =s0.7 rad nT^ (NPCb). The results 
are shown in Fig.[^ respectively panel (a) and panel (b). These 
plots show that LOFAR can provide good constraints for both 
values of the overall extragalactic Faraday depth if a few thou¬ 
sand of lines of sight are used, with better performance for larger 
(Tg. The posterior distributions are much narrower than the prior 
distributions and the posterior means agree with the correct val¬ 
ues within l-2cr. Tests performed with a smaller number of lines 
of sight are presented in Appendix]^ These plots show that the 
algorithm performs well even for Aios ~ 1000 even if, as ex¬ 
pected, the posterior distributions are wider. 


3.2.2. SKA 

The SKA is expected to observe the entire Southern sky with 
a spatial resolution of 2" and a sensitivity in polarization of 
=» 4/rJy/beam (see, e.g., |Johnston-Hollitt et al.||2015) . The re¬ 
sulting sky grid of Faraday rotation values will be 200-300 times 
denser than the largest catalog available at the moment (see e.g. 
Hales |2()T3 i. The better resolution of 2", compared to the 45" 
of Taylor et al. ( 2009| l, will make it possible to identify opti¬ 
cal counterparts uniquely and hence to assign a redshift estimate 
to a larger number of sources through spectroscopic follow-up 
observations. The SKAl Re-Baseline Design 2015 indicates the 
frequency bands 1, 2, and 5 as available on SKA_MID during 
SKA-Phase 1; 

- band 1, 0.35-1.05 GHz; 

- band 2, 0.65-1.67 GHz; 

- band 5, 4.6-13.8 GHz. 

In the following we consider the frequency range 0.65-1.67 GHz, 
since receivers in band 2 (B2) should be constructed first. 

We produce mock Faraday depth values assuming a maxi¬ 
mum standard deviation in the noise distribution of the rota tion 
measure of 0.8 rad m“^, according to Stepanov et al. ( 2008| l, for 
this frequency range and for S /N > 5. We generate a catalog 
of coordinates in the South Polar Cap (SPC) based on the as¬ 
sumption of one polarized source per square degree and we se¬ 
lect those with Galactic latitude b < -55°. This translates in 
Aios ~ 3500. We will refer to this catalog as B2SPC1. We pro¬ 
duce catalogs of Faraday depth values assuming cTg xH.Q rad m“^ 
(B2SPCla) and cTe ;^0.7radm-2 (B2SPClfe). 


rmsf is the half-power width of the Faraday depth spread function. 


In Fig. 1^ we show the results for B2SPCla and B2SPC1/7, 
in panel (a) and (b) respectively. The plots for a-g x!7.0radm“^ 
indicate that with Aios ~3500 we obtain a good disentangling of 
the environmental and intrinsic contributions, and their redshift 
and luminosity dependence, with narrow posterior distributions 
and mean values within sslcr from the true values. For the same 
amount of lines of sight but a smaller overall Faraday rotation 
(Tg xQ.l rad nT^, the dispersion in (Tenv,o, Alum, and Ared becomes 
broader but with mean values still within a few standard devia¬ 
tions from the assumed value. On the other hand, we do not get 
good constraints for the parameter (Tint,o that appears character¬ 
ized by a large dispersion. For a better constraint of this param¬ 
eter we would need to resort to a larger number of lines of sight. 
Since the sensitivity of the SKA will allow to detect hundreds of 
sources per square degree, narrower posterior distribution of the 
parameters can be obtained by increasing the number of lines of 
sight. 

Because of the different number of lines of sight and differ¬ 
ent frequency bands we can not directly compare these results 
with the results for LOFAR. For an analysis of SKA perfor¬ 
mance at lower frequencies see Appendix In this appendix, 
results for SKA observations in band 1 as well as in band 2 but 
with a smaller number of lines of sight are shown. We do not 
consider frequency band 5. Even if moving to higher frequency 
sources have a higher degree of polarization, the maximum un¬ 
certainty in Faraday depth is larger and the total source counts 
reduce. Consequently, we expect that the parameters are poorly 
constrained if the same number of lines of sight is taken. 


3.2.3. ASKAP 


The Polarisation Sky Survey of the Universe’s Magnetism 
(POSSUM, [Gaensler et al. pOlO i is planned between 1130 
and 1430 MHz with the Australian Square Kilometre Array 
Pathfinder (ASKAP). The survey will reach a sensitivity in U 
and Q of <10pJy/beam and a resolution of 10". This will result 
in a density of 70 polarized sources per square degree ( |Hales et| 
al.|2014| l and a Faraday depth catalog one hundred times denser 
than those currently existing (e.g., Taylor et al.|2009| l. To test the 
capabilities of this survey in constraining extragalactic magnetic 
fields, we use the same set-up described in Sect. 3.2.2| but with 
an uncertainty drawn from a Gaussian with CTnoisg = 6.radm“^. 
This should be a reliable approximation of the maximum uncer¬ 
tainty in Faraday depth for a polarized signal with a S fN > 5 
( [Stepanov et al.|2008| l. The results are shown in Fig.[^ ASKAP 
observations appear already quite promising for disentangling 
the contribution intrinsic to the source from that due to the 
large-scale environments for o-g :«7.0radm“^, already if Faraday 
depth values are available only for one source per square degree. 
The dispersion in the parameters is larger than that derived in 
Sect. [322] for the same collection of sources and overall extra¬ 
galactic Faraday variance (Fig.|^), as expected being the max¬ 
imum uncertainty in Faraday depth larger. However, the mean 
values of the posterior agree within at most about Icr with the 
true values. 


3.3. Discussion 


We developed a Bayesian algorithm to disentangle extragalac¬ 
tic contributions in the Faraday depth signal. It builds upon 
an algorithm to reconstruct the Galactic Faraday screen and its 


uncertainty-correlation structure previously presented by Opper- 


jmann et al.j ( |2015j l. We tested the algorithm by modeling the 
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overall extragalactic contribution as the sum of an intrinsic and 
an environmental component, as described in Sect. [m but see 
Appendix [C] for the inclusion of a constant or a latitude depen¬ 
dent term. These tests show that the algorithm is able to discrim¬ 
inate among the different components and their dependence on 
the luminosity and the redshift of the source, if a suitable number 
of lines of sight is available. 

For test purposes, we built mock catalogs according to the 
specifications of the catalogs available in the literature after cor¬ 
recting for poorly known uncertainty information (see Sect. [ 33 . 
The results in Fig.j^b) indicate that by applying the algorithm to 
the currently available dataset we could already infer preliminary 
information about extragalactic magnetic fields, if a few thou¬ 
sands of sources with reliable observational uncertainty were 
available. For the 4003 sources in Fig.j^b), indeed, redshift and 
flux density values are available. Nevertheless, since actually we 
do not consider all their Faraday depth noise variances to be 
reliable (see Appendix A of Oppermann et al. 2015 1 , and the 
applied correction factors, 77,, are only statistical estimates, we 
would need to additionally sample in the 77 ,-parameter space in 
order to use the algorithm with real data , as described by |Op-| 
Ipermann et al.| ( |2015[ ). The modern techniques of analysis, such 
as model fitting of the fractional polarization components q and 
u (e.g., Farnsworth et al.||201 1[ [Ideguchi et al.||2014|l, rotation 
measure synthesis (e.g., Brentjens & de Bruyn|2005 Akahori ^ 
|al.|2014b| ) and Faraday synthesis ( Bell & EnBlin|2012| l, and the 
future radio surveys will partially overcome this problem. The 
large bandwidth of the new radio interferometer will allow to re¬ 
duce the risk of 77;7r-ambiguity, particularly strong when the T^-fit 
approach is used, as well as to reach a sufficiently high resolu¬ 
tion in Faraday depth to distinguish nearby Faraday components. 
We note that when the distance between two peaks in a Faraday 
spectrum is smaller than the resolution, the uncertainty in the 
Faraday depth may be driven by the Faraday point spread func- 
tion (e.g.,|Farnsworth et al.|201 l[|Farnes et al.|2014b[ Kumazaki 


et al. 2014| l. For these reasons, the catalogs coming from Faraday 
depth grids planned with LOFAR, ASKAP, and the SKA will 
be more reliable both in terms of Faraday depth values and of 
uncertainties. Therefore, we investigated the prospects of these 
simpler future datasets here and address the more complex ap¬ 
plication of our technique to present data in a separate work. We 
assumed the computation to be dominated by the uncertainties 
in the Faraday depth estimates for each source, while the co¬ 
ordinates, luminosities and redshifts to be exactly known. With 
the advent of highly-accurate Faraday depth catalogs, this may 
change and in future work we will need to investigate to what 
depths and accuracy redshifts, and luminosities are needed in or¬ 
der to apply the method we propose. This is particularly impor¬ 
tant for redshifts. Indeed, the catalog used in this work reports 
a spectroscopic redshift for each source ( Hammond et al.|2012| l. 
Recently, optical surveys, such as 2MASS, WISE, and Super- 
COSMOS, have measured the photometric redshifts for millions 
of galaxies up to z « 0.5 (e.g. Bilicki et al.|2014) l, and the sur¬ 
veys planned with the next generation of telescopes will further 
increment this number. Nevertheless, photometric redshifts are 
less accurate than spectroscopic redshifts. This will require to 
evaluate the impact of their uncertainties on our results. 

Our tests indicate that LOFAR, ASKAP, and the SKA will al¬ 
low us to infer information about cosmological magnetic fields 
already with a few thousands of lines of sight, with better per¬ 
formance when lower frequencies are used. We want to stress 
that the uncertainties used in our tests only represent statistical 
uncertainties and any systematic issues have been neglected. In 
principle, LOFAR observations could be already used for this 
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aim but the development of a pipeline for the reduction of po¬ 
larization data is still in progress. The enhanced capabilities of 
the Jansky Very Large Array (JVLA) make the new centimeter- 
wavelength sky survey (VLASS, https://science.nrao. 
edu/science/surveys/vlass) planned with this instrument 
a good opportunity for the application of this algorithm and for 
delivering significant results in the study of cosmic magnetism 
already in the next years. 

In Table we give a quantitative summary of the results. 
Lor each test we report the true values r of the © parameters 
which describe the intrinsic and environmental contribution to 
the Laraday depth, their mean values 77 , and their uncertainties 
cr, as well as the displacement of the mean from the true value 
in terms of the uncertainty. 


cr 


These values have been computed after discarding the bum-in 
samples by visual inspection. The comparison of the results ob¬ 
tained with the mock catalogs created for present instraments 
2C1 and 2C2 indicates that by increasing the number of sources 
with known redshift by a factor ten, the uncertainty in the 0 
parameters can be reduced by about a factor two. This would 
not longer be valid if the sources without redshift information 
represent a different population than the observed one, since in 
our test we adopted mock values for these sources randomly ex¬ 
tracted from the observed redshift distribution. A similar result is 
obtained if the observational uncertainty cTaome is reduced by ap¬ 
proximately a factor ten, as shown by the comparison of results 
corresponding to catalog 2C2 and the SKA catalog B2SPCla. 
These two tests refer to a similar number of lines of sight and 
to a similar frequency rang^ In Lig.j^we present a visual sum¬ 
mary of the results for all the two- and three-component sce¬ 
narios corresponding to an overall extragalactic Laraday rotation 
of approximately both 7 rad nT^ (left panels) and 0.7 rad nT^ 
(right panels) for the four parameters (Tint 0 , cTenv.o^ Tium, and;^fred■ 
The uncertainty on the parameters 0 decreases by increasing the 
number of lines of sight, Aios, and by decreasing the observa¬ 
tional error, (T„oise- The best results are obtained when a good 
compromise between these two numbers is used, as for exam¬ 
ple shown by the scenarios NPCa and NPCh. Nevertheless, the 
comparison of the results of tests corresponding to different in¬ 
struments can be not straightforward because they refer to dif¬ 
ferent frequency bands and to slightly different values of the 0 
parameters. 

The problem we are tackling is characterized by different 
complexities. We are looking for a very weak signal, by using 
the residual information left in the data after the Galactic contri¬ 
bution has been derived. At the same time, we expect this signal 
to be characterized by a possible redshift dependence, even if 
very weak since not yet detected, and by cross-correlations of the 
extragalactic magnetic field along different lines of sight. In this 
work we try to address the first two issues with a Bayesian ap¬ 
proach able to properly combine all the available observational 
information and allowing for a redshift dependence in our model. 
0 n one hand, we might overestimate the extragalactic contribu¬ 
tions, since we are not including an uncorrelated Galactic com¬ 
ponent in the model. This term can be easily included in our 
algorithm, as shown by the tests in Appendix and quantified 
with real data. Gn the other hand, we do not take into account 


^ Indeed, the 2C2-sources all belong to the catalog by Taylor et al. 
P009|l, therefore their noise uncertainties refer to the frequency range 


1365-1435 MHz. 
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a possible correlated component of the extragalactic magnetic 
field, which if present would be beneficial for our analysis. This 
means that the magnetic field values that can be inferred have to 
be considered rather as a lower limit, even if we expect this error 
to be small (e.g., Akahori & Ryu|[2011| l. This is not a problem 
for surveys performed with LOFAR, since the expected density 
of extragalactic polarized sources is lower than the density of 
sources we assume. On the contrary, with ASKAP and the SKA, 
we will be able to detect at least one hundred sources per square 
degree. An approach to deal with possible cross-correlations be¬ 
tween lines of sight is being developed with the intention to 
apply this new technique to forthcoming ultra-deep JVLA ob¬ 
servations from the CHILES Con Pol survey (Hales et al. in 
preparation). The combination of the two methods would allow 
to exploit in an optimal way the information coming from the 
denser grids of Faraday depth measurements provided by SKA 
and ASKAP surveys. 

Finally, we stress that the aim of this paper is to give a proof 
of concept. Indeed, the present version of the algorithm requires 
a high computational burden that limits its efficiency. For exam¬ 
ple, the time required to run a test on our computer cluster is of 
the order of a couple of months per chain even when only a few 
thousands lines of sight are considered. We are currently work¬ 
ing on performance optimization to make the algorithm suitable 
to be applied to the upcoming large Faraday depth catalogs. 


3.4. Future developments 

The work presented in this paper prepares a Bayesian tech¬ 
nique to investigate of magnetic fields in the large-scale struc¬ 
ture, in particular in filaments and voids. As a next step we 
want to discriminate among the amount of Faraday rotation due 


to each large scale structure environment (see also Vacca et al. 
2015|l. When different large-scale environments are considered. 


the variance in the extragalactic Faraday rotation can be param¬ 
eterized as 








Li 

Lq] (1+z,)" 


int.O 


^env 

+ 

i=i 




(30) 


where crj, cr 2 ,...,are the contributions from Aenv different 
environments and is the length of the line of sight, i, through 
each environment, j. 

With a Bayesian approach, |Jasche et al.| P010| l (see also 
|Leclercq et al.|[2015[ ) reconstructed the cosmic density field. 
They used optical data from the SDSS Data Release 7 (|^ aza-| 
[jian et aLl|2009| l and classified the structures as voids, sheets, 
filaments, and galaxy clusters, according to the classification 
scheme of |Hahn et al.| ( [2007] l. This reconstruction enables us to 
compute the path covered by the radio signal through each en¬ 


vironment, i.e., the elements /,y in Eq. (30i, once the position of 


the radio source is identified by using the redshift. This posterior 
of the large scale structure density field is available in the form 
of samples. Radio sources can belong to different environments 
(e.g. galaxy clusters and filaments) and the path covered by the 
signal in each environment differs for each radio source. These 
facts can be statistically taken into account with the use of a col¬ 
lection of sources distributed over all the sky and of different 
realizations of the large scale structure. We plan to use the full 
posterior of the large scale structure in order to statistically esti¬ 
mate which is the amount of variance due to the different types 
of environments in the observed Earaday depths. 


4. Conclusions 

The properties of cosmic magnetic fields constitute outstanding 
questions in modern cosmology. To get a better understanding it 
is essential to shed light on the properties of magnetic fields in 
large-scale environments, meaning filaments and voids, where 
turbulent intracluster gas motions have not yet enhanced the 
magnetic field and, consequently, the magnetic field strength and 
structure still depend on the seed field power spectrum. 

Upcoming generations of radio telescopes, first LOEAR, and 
in the next decades the SKA, will perform polarization sky sur¬ 
veys with high sensitivity. Modern techniques based on rotation 
measure synthesis and Faraday synthesis, will enable us to per¬ 
form a proper analysis of the polarization properties of extra¬ 
galactic radio sources, thus providing unprecedented, highly ac¬ 
curate Faraday depth catalogs in frequency ranges from a few 
hundreds MHz to a few GHz. A statistical approach is required 
to exploit the information encoded in these data. For this reason 
we developed a Bayesian algorithm able to combine radio ob¬ 
servations with luminosities and redshifts of sources, aiming at 
disentangling contributions to the extragalactic Faraday rotation 
intrinsic to radio sources and due to the large-scale structure, 
and in this way infer information about large-scale magnetic 
fields. Knowledge of the redshift is essential in this approach. 
The present all-sky photometric optical surveys and the surveys 
planned with the next generation of telescopes will greatly en¬ 
large the number of sources for which this information will be 
available. 

The work described in this paper is a proof of concept and 
shows that our algorithm can be used to discriminate between 
the Faraday depth generated by the radio source itself and the 
contribution due to the large-scale structures. Additionally, our 
algorithm is able to investigate the dependence of these terms 
on the redshift and the radio luminosity of the sources. The tests 
performed with mock LOFAR, ASKAP, and SKA data suggest 
that this technique is promising for the investigation of magnetic 
fields with strengths of a few pQ down to a few nG, when uncer¬ 
tainties in the data are up to a few rad m“^ and known with high 
accuracy. We note that our modeling does not take into account 
any correlated component of extragalactic magnetic fields. Con¬ 
sequently inferred magnetic field strengths have to be considered 
as a lower limit. 

The main characteristics of upcoming polarization surveys 
can be summarized by the number of lines of sights and by the 
maximum observational uncertainty in Faraday depth. Our tests 
indicate that, for a given number of lines of sight, better con¬ 
straints can be obtained with observations at lower frequencies, 
because of the smaller observational uncertainty. Therefore, in 
principle, LOFAR and, thanks to its higher sensitivity even more, 
SKA_LOW (50-350 MHz) observations would be ideal. Never¬ 
theless, the scant number of polarized sources at these frequen¬ 
cies and the difficulties in the calibration of the data could make 
the use of these data complex. ASKAP and SKA_MID obser¬ 
vations respectively in the frequency range 1130-1430 MHz and 
650-1670 MHz appear to be promising as well. We should be 
able to put useful constraints on large scale magnetic fields al¬ 
ready with Faraday depth measurements for a few thousands of 
sources, and improve their determination by increasing the num¬ 
ber of lines of sight. An increment in the number of lines of sight 
by a given factor reduces the uncertainty in the estimation of the 
intrinsic and environmental contribution as a reduction by the 
same factor in the observational uncertainty does, indicating that 
deeper observations of small fields could be a valuable or even 
better alternative to all sky surveys. 
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We are aware that many aspects of our approach require im¬ 
provements: e.g., computational efficiency, inclusion of a corre¬ 
lated extragalactic magnetic held component, and of uncertainty 
in redshift, etc. Nevertheless, we present a hrst step toward a 
Bayesian study of magnetic helds associated with the cosmic 
large-scale structures. 
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‘^int,0 rad/m^ ^env,0 rad/m^ Xhrm Xred 


Fig. 2. Results obtained with a two-component scenario for 41632 (2C1) lines of sight in panel (a) and for 4003 (2C2) in panel (b). In each panel 
the top plots of each column show the 1-dimensional projection of the posterior as well as the true value (dotted line), the outcome of the analysis 
(dashed and dashed-dotted lines), the prior (continuous line). The panels in color show the 2-dimensional marginalized views of the posterior as 
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(b) 




Fig. 3. As Fig.j^but for results obtained with a two-component scenario for LOFAR HBA observations for Nias ~2200 and an overall extragalactic 
standard deviation in Faraday depth of (a) x: 7 radm^^ (NPCa) and (b) 0.7 radm^^ (NPCZ?). 
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Fig. 4. As Fig.j^but for results obtained with a two-component scenario for SKA observations in the frequency range 650-1670 MHz for Aios ~3500 
and an overall extragalactic Faraday rotation of (a) « 7.0radm“^ (B2SPCla) and (b) « 0.7 radm^^ (B2SPC17>). 
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Fig. 5. As Fig.j^but for results obtained with a two-component scenario for ASKAP observations in the frequency range 1 130-1430MHz for 
Aios ~3500 and an overall Faraday depth of «7.0radm“^ (POSSUM). 
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Fig. 6. Uncertainty on crint,o, crenv,o,;t'ium, and;^red (from the top to the bottom) as a function of the observational uncertainty (Xnoise and the number 
of lines of sight for an overall extragalactic Faraday depth variance of «7.0radm“^ (left panels) and radm“^ (right panels). The size 
of the points is proportional to the uncertainty on the parameter. A description of this uncertainty is given also by the greyscale. Circles refer to 
scenarios including two components, while triangles refer to three-component scenarios. Continuous lines correspond to <^noise/^ios =const. 
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Appendix A: Extragalactic Faraday depth variance 

In this appendix we derive Eq. Q from Eq. (5). In Eq. Q we define the extragalactic Earaday depth variance as 




r'' di n di' 
Jo dzJo dz' 


{n,(z)n,(z')B,(z)Bi(z')}dzdz'. 


(A.l) 


This definition can be expressed also as a function of distance along the line of sight I 

-Kzi) H(Zi) 




riKZi) 

di 

Jo Jo 


dl'{n,(DnS')Bi(l)Bi{l')). 


(A.2) 


If we assume that the thermal gas density is not characterized by significant fluctuations and define a new variable r - I' - I, we 
have 


{nMne(l')Bi(l)Bi(l')) ^ {nl{l)){B,{l)B,{l + r))(BM.) = («e(0)Cij(r|ne), 


(A.3) 


where CB(r|ne) is the conditional magnetic field correlation function for an environment with thermal gas density n^. Indeed, we 
expect the magnetic field strength to be a function of the thermal gas density. With these new definitions, the extragalactic Earaday 
depth variance reads 
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dr(n^Jl))CB(rln,), 


which can be further simplified if we consider the limit of a infinitely far away source. 
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Recalling the definition of correlation length, 
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(A.5) 


(A.6) 


(A.7) 


As described in Sec.|^ in a homogeneous Universe tig - neo(l + z)^, Ai = Aio(l + z) *, and {B^} = (Rq)(1 + z)*. Therefore, 

dl 




j-4 r 

Jo 


<ng>Aoi<<) (l+z)5dz, 
dz 


(A.8) 


where the increment in wavelength due to the expansion of the Universe has been taken into account as well and we have assumed 
z(l) ~ z(l + r). If we use the definition of proper displacement along a light-ray d//dz = c(l -H z)~^ IH{z), this finally leads to Eq. (j^. 


( 0 ; 


u =4 r 

Jo 


(n^^}Ao,(B^o,}—(l +zM. 

In an inhomogeneous Universe with different environments, the differential variance in each environment is 


(A.9) 


= {nl{z)){B\B\nA 

dx 

Therefore, it follows 

dx dz 


( 0 , 


, , r^' dv dz r c 




RM 

dx 


(A. 10) 


(A. 11) 


where the integral can be replaced by a discrete sum over typical environments (see also Eq. (|30ll). 


Article number, page 17 of|27| 








A&A proofs: manuscript no. paper 


Appendix B: Convergence 

We start the algorithm from two random positions in the 0-space and then explore the space ^ until convergence. Each of these 
two sequences of steps in the 0-space is called Gibbs chain. To assess convergence of each Gibbs chain, we require the following 
conditions to be satisfied: 


- the number of steps taken from each chain to be at least about 101^ for each parameter, where Ic is the number of steps at which 
the correlation coefficient drops to 10%; 

- the Gelman and Rubin test ( jGelman & Rubin |1992|[Brooks & Gelman |1997| l. We evaluate the intra-chain variance 

(B.l) 

and the inter-chain variance 

(B.2) 

m - 1 ^ 

for our 0 parameters and use them to compute the potential scale reduction factor R defined as 



where m is the number of chains, n is the half-length of each chain, 6j and Sj are respectively the mean and the standard deviation 
of the y'th chain, while 6j is the mean of the chain j and 9 is the variance of the chain means. We require R - \ within a few 
percents for each parameter. 

We consider the chains to be converged when both these conditions are satisfied. 



Fig. B.l. Potential scale reduction factor R evaluated every 10 samples for the parameters (Tint.o, cTem.o, Xivim, ^.ndxrei, as a function of the number 
of steps in the MCMC. The continuous and the dashed lines represent a potential scale reduction factor equal to 1.0 and 1.1, respectively. It can be 
seen that the Gelman and Rubin test indicates convergence after typically a few thousand steps. 


As an example, in Fig. |B.l| we show the plot of the potential scale reduction factor R versus the number of steps for the 0 
parameters in the scenario 2C1. The potential scale reduction factor has been evaluated every 10 samples. After about 5000 steps 
we obtain R = 1 within a few percent. In Table B.l the potential scale reduction factor R for each parameter of each test is reported. 
For present-instrument (scenario 2C1), FOFAR (scenario NPCa) and SKA (scenario B2SPCla) observations, in Fig. |B.2| we show 
the correlation coefficient pj 


Pjik) = 


Ztl'jXs - 0j)(Xsrk - Oj) 


(B.4) 


as a function of the number of steps for each parameter (Xenv.o^ Tium, Here Xs is the value of one of these parameters 

for a given step. For all the parameters in the scenario 2C1, /c turns out to be ^3000. while for the scenarios NPCa and B2SPCla, 
the correlation length is about 500-1000. In these plots the steps in the burn-in phase have been discarded by visual inspection. 


^ We tuned the variance of the Gaussian step proposal to ensure an acceptance rate of the Metropolis-Hustings algorithm of approximately 15-30% 
Article number, page 18 of 1 


27 




















V. Vacca et al.: Using rotation measure grids to detect cosmological magnetic fields - a Bayesian approach 





Fig. B.2. Correlation coefficient p as a function of the number of steps in the MCMC for the four © parameters in the 2C1 test (top and middle 
panels) and for the crint,o parameter in the LOFAR NPC and the SKA B2SPClfl tests (bottom panels). The steps in the burn-in phase have been 
discarded. The dashed line indicates a correlation coefficient of 0.1, while the continuous line the zero-level. 


Appendix C: Alternative scenarios 


In this appendix we present the application to scenarios including three components, representing: 


- an intrinsic, an environmental, and a constant contribution (scenario 3C), 




Uoj (l+z,)4 


+ 


^i(^hXred) ^2 

a; 


+ cr 


2 

c 


(C.i) 


The constant contribution cr^ takes into account terms that are not described by the parameterization of the other two but that 
could nevertheless be present in our data (e.g. the ionosphere, under the assumption that this does not show any direction- 
dependence); 

- an intrinsic, an environmental, and a latitude-dependent contribution (scenario LD), 




T XA'Ium 

.0) = !^ 


int,0 


LoJ (l+z,)4 


, Di(Zi,Xce6) 2 , n \ Z 

+ --CTenv.O + T’(^)c^lat 


(C.2) 


where p{b) is the Galactic profile from Oppermann et al. (2015|. The latitude-dependent contribution p{b)(T^^^ may explain a 
residual latitude dependence not taken into account in the modeling of Oppermann et al. P015| l, e.g. an uncorrelated Galactic 
signal not captured by their analysis. 
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Table B.l. Potential scale reduction factor R from the Gelman and Rubin test. 


ID 

^int,0 

^ env,0 

zVc/lat 

/^lum 

Xred 

2C1 

1.00 

1.01 


1.00 

1.01 

2C2 

1.00 

1.01 


1.00 

1.01 

3C1 

1.01 

1.12 

1.07 

1.00 

1.09 

3C2 

1.00 

1.00 

1.00 

1.01 

1.00 

LDl 

1.03 

1.05 

1.00 

1.02 

1.04 

LD2 

1.00 

1.01 

1.01 

1.00 

1.01 

GW 

1.00 

1.00 


1.00 

1.00 

NPCa 

1.00 

1.00 


1.00 

1.00 

POSSUM 

1.02 

1.03 


1.01 

1.03 

B2SPClfl 

1.00 

1.01 


1.00 

1.01 

B2SPC2 

1.00 

1.00 


1.00 

1.00 

BlSPClfl 

1.00 

1.00 


1.00 

1.00 

B1SPC2 

1.00 

1.00 


1.00 

1.00 

PO 

1.02 

1.03 


1.01 

1.03 

PI 

1.01 

1.02 


1.00 

1.02 

NPCi^ 

1.00 

1.00 


1.00 

1.00 

B2SPCli 

1.45 

1.08 


1.01 

1.03 

BlSPClfo 

1.04 

1.02 


1.01 

1.01 


As for the simplest 2C-scenario, for each of these scenarios we run two tests corresponding to 41632 and 4003 lines of sight. Fig. |C.l| 
and Fig. C.2 show the results respectively when a third constant component (tests 3C1 and 3C2) and latitude-dependent component 
(tests LDl and LD2) are included. In both figures, we show the plots for 41632 lines of sight in panel (a) and for 4003 lines of sight 
in panel (b). 

These plots indicate that the algorithm performs well also when three components are considered. As expected, the values of 
the © parameters recovered are less accurate when a lower number of lines of sight is used. The additional parameters tend to lead 
to a slight increase in the posterior uncertainty for the other parameters when comparing with the results of 2C1 and 2C2. 


Appendix D: Priors 

In order to get a data-driven solution and to keep our assumptions as general as possible, an uninformative prior should be adopted 
for the 0 parameters. Since we included in our model all the main redshift and luminosity dependencies, the Gaussian prior in 
Eq. (j2T|i is suitable for;^flu^l and;^^red- Concerning cri„( 0 and we may ask if different priors can have an impact on our results. 

In Sect. |3.1[ we adopted a flat-prior, 

P(cr) - const. (D.l) 

In this appendix we present two tests corresponding to extreme choices of these priors. Indeed, we considered a flat prior in cr^ 
(scenario PI) 

P(cr^) - const, (D.2) 


and a flat prior in In(cr^) (scenario PO) 
P(ln(cr^)) = const. 


(D.3) 


In Fig. 


D.l 


we show the results for the two priors; flat in cr^ in panel (a) and flat in ln((T^) in panel (b). 


The first choice is an optimistic prior, since it implies a suppression of cr values << 1, pushing for the recovery of larger, 
possibly ~ 1 values of cr. The second choice is a pessimistic prior since it would weight all small and large cr values in the same 
way, favoring negligible extragalactic contributions, easily compatible with the data due to the shape of the likelihood. We stress 
that the final results are not affected by the choice of the prior. However, this choice has an impact on the convergence-time, since 
the starting point of each chain is randomly extracted from the prior. For example, for a flat prior in In(cr^), the chances to extract a 
very small cr value (<< 1) are larger than for the prior used in Sect. [m possibly making the convergence time very long. For theses 
tests, the number of lines of sights and the assumed noise properties are the same as in the 2C1 test. 


Appendix E: LOFAR 


An interesting region of the sky is represented by the Great Wall, where one of the largest filaments of optical galaxies has been 
observed ( |Gott et aL]|2005| l. We generate a mock collection of source coordinates in a region of the sky as large as this region 
(7.5*’ <RA< 17.5*’ and 25° <Dec< 65°), considering a density of one polarized sources per 1.7 square degrees (survey of 8h 
per pointing) and Galactic latitude b > 55°. This results in Aios ~ 1000 and a maximum uncertainty in Faraday depth CTnoise = 
0.05 rad m“^. We will refer to this moc k catalog in the following as GW. The results for an overall extragalactic Faraday rotation 
of ^ 7.0 rad m“^ are shown in Fig. E.l and indicate that it is possible to disentangle the intrinsic and environmental contributions 


already with a number of lines of sight of about one thousand. 
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Appendix F: SKA 

In this appendix we investigate the performance in the frequency band 1 (Bl), meaning 350-1050 MHz, with the same collection 
of sources (approximately 3500) described in Sect. |3.2.2| According to |Stepanov et al.| ( |2008[ ), this frequency range corresponds to 
an expected maximum uncertainty in Faraday depth of 0.3 rad m“^, for a S/N > 5. We will refer to this catalog as BlSPCl. We 
generate mock rotation measure catalogs con sidering an overall extragalactic Faraday rotation of ^ 7.0rad m“^ (BlSPCla) as well 


as (Te ~ 0.7radm (BlSPClh). In Fig.F.l we show the results. These tests indicate that a catalog in band 1 allows derivation 


of better constraints then a catalog in band 2. This becomes particularly evident for an overall extragalactic Faraday rotation of 
ss 0.7 radm^^. 

Additionally, we examine the case when a smaller number of lines of sight is available. Therefore, we generate a catalog 
corresponding to a density of one polarized source per three square degrees, including all the sources with Galactic latitude b < -55°. 
This translates in Aios ~ 1000. We will refer to these catalogs as B1SPC2 and B2SPC2, respectively for the frequency ranges 
corresponding to band 1 and band 2. We generate these catalogs assuming an overall extragalactic Faraday rotation of ss7 radm^^. 
In Fig. F.2 we show the results for both frequency bands. 

When data in the frequency range from 350 to 1670 MHz are used, already a number of lines of sight of the order of one 
thousand is enough to disentangle the contribution intrinsic to the source from the one due to the environment between the source 
and the observer for cTe 7.0rad/m^. Our results indicate that the posterior distributions narrow moving to lower frequencies with 
the best performance obtained in the frequency band 350-1050 MHz. At higher frequencies the posteriors are broader but with mean 
values still in agreement with the real values. 
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(b) 




*^int,0 rad/m^ *^env,0 rad/m^ *^lat rad/m^ Xlum Xred 


Fig. C.l. Results obtained with a three-component scenario, including a constant contribution, for 41632 (3C1) lines of sight (a) and 4003 (3C2) 
lines of sight (b). 
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*^int,0 rad/m^ *^env,0 rad/m^ *^lat rad/m^ Xlum Xred 

Fig. C.2. Results obtained with a three-component scenario, including a latitude-dependent contribution, for 41632 (LDl) lines of sight (a) and 
4003 (LD2) lines of sight (b). 
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(b) 




Fig. D.l. Results obtained with a two-component scenario and an overall extragalactic Faraday 7 radm ^ for (a) a flat prior in cr- (PI) and for (b) 
a flat prior in ln((T^) (PO) for 41632 lines of sight. 
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Fig. E.l. Results obtained with a two-component scenario for LOFAR HBA observations for an overall Faraday depth of »:7.0radm ^ and 
^1000 (GW). 
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Fig. F.l. Results obtained with a two-component scenario for SKA observations in the frequency range 350-900 MHz for Nias ~3500 and an 
overall Faraday depth of (a) ^T-Oradm^^ (BlSPCla) and (b) ^O-Tradm^^ (BISPCIZ7). 
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Fig. F.2. Results obtained with a two-component scenario for SKA observations in the frequency range (a) 350-900 MHz (B1SPC2) and (b) 
650-1670MHz (B2SPC2) for Nios ~1000 and an overall Faraday depth of sjT.Oradm^^. 
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